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Abstract 

Exciton effects on conjugated polymers are investigated in soliton lattice states. 
We use the Su-Schrieffer-Heeger model with long-range Coulomb interactions. 
The Hartree-Fock (HF) approximation and the single-excitation configuration- 
interaction (single-CI) method are used to obtain optical absorption spectra. 
The third-harmonic generation (THG) at off-resonant frequencies is calculated 
as functions of the soliton concentration and the chain length of the polymer. 
The magnitude of the THG at the 10 percent doping increases by the factor 
about 10 2 from that of the neutral system. This is owing to the accumulation 
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of the oscillator strengths at the lowest exciton with increasing the soliton 
concentration. The increase by the order two is common for several choices of 
Coulomb interaction strengths. 
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1. Introduction 



Recently, effects of electron-electron interactions have been investigated in 
many aspects of conjugated polymers. For example, electronic excitation struc- 
tures in the half-filled conjugated polymers with the constant bond alternation 
have been theoretically investigated by using the exciton formalism [1] and 
the exact diagonalization method [2], and also by solving the time-dependent 
Hartree-Fock (HF) equations [3]. The roles of excitons were pointed out ear- 
lier, but excitation structures have been considered intensively still recently, 
relating with origins of nonlinear optical spectra [2-4]. The lowest energy exci- 
tation has the largest oscillator strength as the most remarkable consequence of 
correlation effects. This is clearly seen when the optical spectra calculated by 
using the HF wavefunctions are compared with the spectra with the correlation 
effects. This feature is observed when the higher correlations are taken into 
account by the single-excitation configuration-interaction (single-CI) method 
[1] and also by the time-dependent HF formalism [3]. 

It is widely known that the soliton, polaron, and bipolaron lattices are 
present [5], when the SSH model [6], its continuum version [7], and the ex- 
tended model with the term of the nondegeneracy [8] are doped with electrons 
or holes. The new bands related with the nonlinear excitations develop in the 
Peierls gap as the doping proceeds. When correlation effects are considered 
by the single-CI, the excitation structures exhibit the presence of excitons. In 
the previous paper [9] , we have considered exciton effects in the soliton lattice 
states of the doped conjugated polymers. There is one kind of exciton in the 
half-filled system, where the excited electron (hole) sits at the bottom of the 
conduction band (top of the valence band). We have called this exciton as the 
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"intercontinuum exciton" [9]. In the soliton lattice states of the doped SSH 
model for degenerate conjugated polymers, there are small gaps between the 
soliton band and the continuum states, i.e., valence and conduction bands. 
Therefore, the number of the kind of excitons increases, and their presence is 
reflected in structures of the optical spectra. A new exciton, which we have 
named the "soliton-continuum exciton" [9], appears when the electron-hole 
excitation is considered between the soliton band and one of the continuum 
bands. We have looked at variations of relative oscillator strengths of the 
new excitons, the soliton-continuum and intercontinuum excitons. When the 
excess-electron concentration is small, the ratio of the oscillator strengths of 
the soliton-continuum exciton increases almost linearly with respect to the con- 
centration. The oscillator strengths accumulate rapidly at this exciton as the 
concentration increases. The contribution from the soliton-continuum exciton 
is more than 90 percent at the 10% doping. 

The purpose of this paper is to investigate how the above changes of the 
characters of optical excitations are reflected in the nonlinear optical proper- 
ties. We might be able to realize large optical nonlinealities in soliton lattice 
systems, because the energy gap is small in doped conjugated polymers and the 
oscillator strengths accumulate rapidly at the soliton-continuum exciton as the 
doping proceeds. Even though experimental facts of the doped polymers have 
not been reported so much, it would be quite interesting to demonstrate the- 
oretically how large optical nonlinearities would be obtained when conjugated 
polymers are doped with electrons or holes up to as much as 10 percent. The 
SSH model with the long range Coulomb interactions of the Ohno expression 
[10] is solved with the HF approximation, and the excitation wavefunctions of 
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electron-hole pairs are calculated by the single-CI. We consider the off-resonant 
nonlinear susceptibility as a guideline of the magnitude of the nonlinearity. In 
this case, multi excitations, such as double (triple) excitations [11], do not 
contribute significantly, so we can assume that the single-CI well describes 
the characters of excitations at off-resonances. We will calculate the third 
harmonic generation (THG) Xtug^) = x(3^;^,^,^) at the zero frequency 
iv — 0, with changing the chain length and the soliton concentration. We will 
show that the magnitude of the THG at the 10 percent doping increases by 
the factor about 10 2 from that of the neutral system. This is owing to the 
accumulation of the oscillator strengths at the lowest exciton with increasing 
the soliton concentration. The increase by the order two is common for several 
choices of Coulomb interaction strengths. 

This paper is composed as follows. In §2, the model is introduced and the 
numerical method is explained. Results of the optical nonlinearity are reported 
and discussion is made in §3. The paper is summarized in the final section. 



2. 



Model 



We use the SSH hamiltonian [6] with the Coulomb interactions: 



H = #ssh + H int . 



(1) 



The first term of eq. (1) is: 



#ssh = -^2(t -ayi)(c} 7<r Ci +lta + h.c.) 




(2) 
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where t is the hopping integral of the system without the bond alternation; a 
is the electron-phonon coupling constant which changes the hopping integral 
linearly with respect to the bond variable yi, c^ c is an annihilation operator 
of the 7r-electron at the site % with spin a; the sum is taken over all the lattice 
sites of the periodic chain; and the last term with the spring constant K is 
the harmonic energy of the classical spring simulating the cr-bond effects. The 
second term of eq. (1) is the long-range Coulomb interaction in the form of 
the Ohno potential [10]: 



flu* = uYSA,^ - ^ L )( c u c *.i ~ ir) 



2 

+ W ( r i,j) E 4,a c i,a ~ n e i) Q2 c] )T c jtT - n el ) , (3) 

where n e \ is the number of 7r-electrons per site, r\j is the distance between the 
ith and jth sites, and 

W(r) = . 1 (4) 

+ {r/aVf 

is the Ohno potential. The quantity W(0) = U is the strength of the onsite 
interaction, V means the strength of the long range part, and a is the mean 
bond length. 

The model is treated by the HF approximation and the single-CI for the 
Coulomb potential. The adiabatic approximation is applied to the bond vari- 
ables. The HF order parameters and bond variables are determined selfconsis- 
tently using the iteration method [12]. The details of the formalism have been 
explained in the previous paper [9]. A geometry of a ring with the coordinate 
of jth carbon atom (1 < j ' < N), 



is used for a polymer chain, in order to remove edge effects. Here, R = 
Na/(2n) is the radius of the polymer ring, N is the system size, and a is the 
lattice constant. The electric field of light is parallel to the x-y plane. The 
optical absorption spectra where light is along with the x- and y-directions are 
summed as shown by the following formula: 

J2 E k P(uj - E K )((g\x\K)( K \x\g) + (g\y\K)( K \y\g)). (6) 

K 

Here, P(ou) = 7/[7t(c<j 2 + 7 2 )] is the Lorentzian distribution (7 is the width), E K 
is the electron-hole excitation energy, \k) is the «th excitation, and \g) means 
the ground state. 

The THG is calculated with the conventional formula [13-15]: 

X {3) {3u;uj,uj,uj) = e 4 N d Y,fg,vfv,\f\,JK, g 

1 



X 



(E^ g -3u;)(E Xig -2u;)(E K , g -u;) 
1 

+ {El g + u){E^ g -2uj){E^-uj) 
1 

+ (E; tg + u)(El g + 2u J )(E Ktg -Lj) 

+ (E;, g + uj)(El g + 2uj)(El g + 3uj) ] > (7) 
where, iVd is the number density of the crystalline polymer, and this is material- 
dependent. For demonstration of the magnitude of the THG, we use the value 
of the number density of the CH unit, which is taken from trans-polyacetylene: 
iVd = 5.24 x 10 22 cm~ 3 [16]. We also use t = 1.8eV in order to look at numerical 
data in the esu unit. In eq. (7), f\ tK = (\\x\k) is the dipole matrix element 
with the electric field parallel with the x axis, E K g = E K — E g , E K is the energy 
of the excited state, and E g is the energy of the ground state. In the actual 
calculation, we change the order of terms so as to take into account the mu- 



tual cancellation among them [17]. Also, we include a small imaginary part in 
the denominator: for example, E Kg — > E Kg + ii] and E* g — > E* — irj. This 
assumes a lifetime broadening, and suppresses the height of the delta-function 
peaks. The THG at uj = does not sensitively depend on the choice of rj. 
This can be checked by varying the broadening. In the next section, we report 
the results with the value rj = 0.021 

The system size is chosen as N = 80, 100, 120 when the electron number is 
even (it is varied from N cl = N, N + 2, N + A, N + 6, N + 8, N + 10 to iV + 12), 
because the size around 100 is known to give well the energy gap value of the 
infinite chain. More larger system size becomes tedious for doing single-CI 
calculations which call for huge computer memories. In principle, we have to 
adjust parameters and find appropriate ones in order to reproduce experimen- 
tal data, such as, the energy gap and the bond alternation amplitude. But, 
we will change parameters arbitrary in a reasonable range in order to look at 
general properties of the optical nonlinearity. We take two combinations of 
the Coulomb parameters (U, V) = (2t, It) and (At, 2t) as the representative 
cases. The other parameters, t = 1.8eV, K = 21eV/A 2 , and a = 4.1eV/A, 
are fixed in view of the general interests of this paper. All the quantities of 
energy dimension are shown in the units of t. 

3. Numerical Results and Discussion 

Figure 1(a) shows the optical absorption spectrum at the 2% soliton concen- 
tration for (U, V) = (At, 2t). The broadening 7 = 0.05t is used. The same data 
have been used in Fig. 3(b) of the ref. [9]. There are two main features around 
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the energies 0.7t and I At. The former originates from the soliton-continuum 
exciton, and the latter is from the intercontinuum exciton. 

Figure 1(b) displays the absolute value of the THG against the excitation 
energy uj. The abscissa is scaled by the factor 3 so that the features in the 
THG locate at the similar points in the abscissa of Fig. 1(a). The small feature 
at about uj = 0.22t comes from the lowest excitation of the soliton-continuum 
exciton and the larger features at about uj = 0.24t and 0.32t come from the 
higher excitations. The features from the intercontinuum exciton extend from 
uj = 0.48t to the higher energies. In the present calculations, the THG in 
the energy region higher than 0.5t is not large relatively. The point, that the 
THG becomes smaller as the excitation energy increases, has been seen in the 
calculations of the half-filled system [4]. However, in the time-dependent HF 
formalism, the THG in higher energies is still larger as shown by Fig. 4 of the 
ref. [3]. The difference of the distribution of the THG strengths might come 
from the difference of the approximation method for electron correlations. 

The THG data like in Fig. 1(b) are calculated for the three system sizes, 
iV = 80, 100, 120, and for the soliton concentrations up to 10%. As clearly 
seen for example in Fig. 1(b), the off-resonant THG at uj = is quite far from 
features coming from excitons. The contributions from double (triple and so 
on) excitations would be very small. Thus, the single-CI calculations could be 
used as a measure of the optical nonlinearities of doped conjugated polymers. 

Figures 2(a) and (b) display the variations of the absolute value of Xthg(0) 
for ([/, V) = (2t, It) and (4t,2t), respectively. The plots are the numerical 
data, and the dashed lines are the guide for eyes, showing the overall behavior 
for each system size. The deviations of the plots from the expected smooth 
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behavior might come from the quantum effect due to the finite system size 
[18]. The linear absorption has the size consistensy, so the plots of the relative 
oscillator strengths vary smoothly as functions of the soliton concentration [9]. 
However, the THG is not size consistent, and spectral shapes depend on the 
system size when N is as large as 100 [18]. Therefore, it would not be strange 
even if |XthgI(0) i s sensitive to the system size and the soliton concentration. 
The THG increases as the system size increases. This behavior is the same as 
has been seen in the calculations of the half-filled systems [18]. 

The increase of the off-resonant THG near zero concentration is very rapid, 
but the THG is still increasing for a few percent to 10% soliton concentration. 
The behavior is similar for the two sets of the Coulomb interactions for Figs. 
2(a) and (b). The THG in Fig. 2(b) is about ten times smaller than that of Fig. 
2(a). This is due to the difference of the interaction strengths. The increase 
by about the factor 10 2 at the 10% doping from that of the neutral system 
seems to be a remarkable fact. This fact might be independent of Coulomb 
interactions strengths characteristic to optical excitation of the system. The 
experiments of nonlinear optics on doped conjugated polymers have not been 
reported so many times. However, the present theory could used as one of 
guidelines for strengths of optical nonlinearities in doped systems. 

Then, why such the large increase of the THG would occur upon doping 
of the polymers? In the soliton lattice theory with the continuum model [5], 
the energy gap decreases as the soliton concentration increases. Therefore, 
it may seem first that the decrease of the energy gap is one of the reasons. 
But, as shown in the Fig. 4 of the ref. [9], the lowest optical gap due to the 
single-CI calculations is almost independent of the concentration, and thus the 
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change of the optical gap would not be the main reason. However, as we have 
discussed in [9], the ratio of the oscillator strength of the soliton-continuum 
exciton increases very rapidly. In view of this change of the exciton characters, 
it would be natural to conclude that the increase of the THG by the factor 
10 2 is due to the fact that the oscillator strengths accumulate rapidly at the 
lowest exciton with increasing the soliton concentration. 

4. Summary 

We have considered the off-resonant nonlinear susceptibility as a guideline 
of the strength of the nonlinearity in the doped conjugated polymers. We 
have calculated the off-resonant THG with changing the chain length and the 
soliton concentration. We have shown that the magnitude of the THG at the 
10 percent doping increases by the factor about 10 2 from that of the neutral 
system. The increase by the order two is common for the several choices of 
Coulomb interaction strengths. 
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Figure Captions 



Fig. 1. (a) The optical absorption spectrum and (b) the absolute value of 
the THG, for the system size N = 100, the electron number iV el = 102, and 
(U,V) = (At,2t). The broadening 7 = 0.05t is used in (a), and rj = 0.02t is 
used in (b). 

Fig. 2. The absolute value of the THG at uo = v.s. the soliton concentration 
for (a) (U, V) = (2t, It) and (b) (4t,2t). The numerical data are shown by the 
triangles (N = 80), circles (N = 100), and squares (N = 120), respectively. 
The dashed lines are the guide for eyes. 
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